Theory of light amphfication in active fishnet metamaterials 
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We establish a theory that traces hght amphfication in an active double-fishnet metamaterial back 
to its microscopic origins. Based on ab initio calculations of the light /plasmon fields we extract 
energy rates and conversion eflaciencies associated with gain/loss channels directly from Poynting's 
theorem. We find that for the negative refactive index mode both radiative loss and gain outweigh 
resistive loss by more than a factor of two, opening a broad window of steady- state amplification 
(free of instabilities) accessible even when a gain reduction close to the metal is taken into account. 
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Recent advances in metamaterials research and active 
plasmonics bring with them the prospect for next gen- 
eration, loss- free optical metamaterial designs [1]. The 
double-fishnet metamaterial, first demonstrated in the 
near infrared [2], is a simple nanoplasmonic structure 
with a broadband effective negative refractive index, the 
origins of which are now well established [3]. To date it 
remains the prime example of a metamaterial able to op- 
erate at visible wavelengths [4] and serves as a template 
for novel optical metamaterial designs [5]. In plasmonics 
an understanding of how losses experienced by surface- 
plasmon polaritons (SPPs) and localized plasmons can 
be overcome by gain media placed adjacent to metals [6- 
8] has had far reaching impact. This is best illustrated by 
the realization of SPP amplifiers [9] and nano-plasmonic 
lasers [10] but also by recent developments in optical 
metamaterials research. In 2010, first evidence has been 
given that loss compensation in active plasmonic meta- 
materials is possible: Experimental [11] and theoretical 
[12, 13] absorption/transmission/reflection (ATR) stud- 
ies on optically pumped gain-enhanced double-fishnet 
metamaterials showed that gain can overcome the inter- 
nal dissipative losses while retaining and in fact enhanc- 
ing the effective negative refractive index response. Fur- 
ther theoretical work has shown that loss compensation 
is achieved when the gap-SPP becomes undamped [14] 
causing the resonance linewidth to be limited by scat- 
tering losses only [14, 15]. The physical origin of light 
amplification in active fishnet metamaterials lies in the 
complex interaction of the strong local fields of the gap- 
SPP mode with the spatially inhomogeneous inversion 
generated by a pump process [12, 16]. Ultimately, to 
develop loss-management strategies for future practical 
metamaterials, a theory is required that traces the in- 
ternal energy flow, inaccessible by experiment, on the 
nanometer scale. 

In this Letter, we present a theoretical framework that 
allows us to describe the microscopic dynamical pro- 
cesses, which, collectively, make light ampliflcation in ac- 
tive metamaterials possible. We conduct, under realistic 
assumptions, extensive numerical pump and probe stud- 
ies and extract together with the macroscopic spectral 
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Figure 1. (color online) Schematic of an active plasmonic 
metamaterial illuminated by an incident plane wave of inten- 
sity / that creates a power fiux into the volume V through 
the surface Az. According to Poynting's theorem [see Eq. (1)] 
the electromagnetic energy U inside the volume changes due 
to energy transfer into and out of various loss (blue) and gain 
channels (red). Note that for normal incidence onto an infi- 
nite film, the fiux between neighbouring volumes cancels. 



response the effective rates of energy exchange between 
physical subsystems of an active fishnet metamaterial. 
The results demonstrate the intricate mechanisms that 
make steady-state amplification in nano-plasmonic struc- 
tures possible: the efficiency of the pump scheme affect- 
ing the spatial deposition of inversion, the electric field 
enhancement that increases the effective gain, and the 
dominance of radiative coupling over internal losses. 

Our studies are based on the Maxwell-Bloch (MB) ap- 
proach [15, 17] combined with a suitable formulation of 
Poynting's theorem. We solve Maxwell's field equations 
in full- vectorial three-dimensional form and in interaction 
with the dynamic polarization responses ascribed to the 
metal and four-level gain medium microscopically, at ev- 
ery point in space, while simultaneously tracing the flow 
of electromagnetic energy into and out of microscopic 
channels (see Fig. 1). This generic method allows us to 
dynamically record and analyze the gain and loss rates 
associated with the physical subsystems of any struc- 
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ture and thereby to identify transient and steady-state 
regimes of operation as well as to precisely determine the 
conditions for amplification and lasing. 

We start by recalling Poynting's theorem, according to 
which the rate of change of the electromagnetic energy 
inside a volume V that contains metal (loss) and four- 
level gain media is determined by 

{du/dt) = -(V • S) - (Pf • E) - (Pa • E) - (Pe • E) (1) 

where u{T,t) = ^ [£o^(r)£^^(r, t) + /io^^(r, t)] is the en- 
ergy density, S(r, t) the Poynting vector, E(r, t) the elec- 
tric field, and Pi(r, t) the polarizations of the free elec- 
tron plasma in the metal {i = f) and the gain media 
at the transitions of absorption and emission {i = a, e). 
The operator (•) performs an integration over V and a 
time- aver aging over At ^ 27r/cj eliminating fast phase 
oscillations. Equation (1) expresses the change of elec- 
tromagnetic energy in V as influx through the surface 
minus the work that the field exercises on the electrons 
in the metal and on the transitions of the gain media. 
It will prove useful to write Eq. (1) as a rate equation: 

dU/dt = -TtU = -AU - r^u - r^U - TeU (2) 

where U = {u -\- ^^Wi) is the total field energy in 
the volume. Ft the energy decay rate, A = (V ■ S) /U 
the flux rate that measures the outflow of energy, and 
Pi = [(Pi ■'E — dwi/dt)]/U the power loss/gain rates. Wi 
are contributions of the dynamic responses to the overall 
electromagnetic energy in the medium that can be calcu- 
lated as in [18]. Note that the so-defined effective rates 
are time-dependent quantities and become constant when 
U is either constant or decreases exponentially during sin- 
gle mode decay. We also emphasize that this formulation 
of volume-averaged energy transfer rates is exact (since 
it directly reflects energy conservation) and applicable to 
arbitrary linear and nonlinear systems. 

The object of our studies is a 260 nm thick fishnet 
metamaterial consisting of two 40nm-thick silver films, 
separated by 60 nm and embedded in an active dielec- 
tric host of refractive index n = 1.62 (see Fig. 1). The 
silver films are perforated with rectangular holes of di- 
mensions ax = 120 nm and ay = 80 nm with square pe- 
riodicity p = 280 nm. For the gain medium (dye) we 
assume a density of = 3.89 x 10^^ cm~^, wavelength of 
maximum absorption (emission) A = 680 nm (710 nm), 
and associated cross-section of a = 7.45 x 10~^^cm^ 
(5.78 X 10-^^ cm^). To model the response of silver we 
use a three-pole Drude-Lorentz model [19]. 

In our analysis, we first pump the active structure with 
a pulse of Ap = 680 nm and then probe it with a signal 
pulse of center- wavelength Ag = 710 nm. Both the pump 
and the probe fields are x-polarized, and the probe falls 
within the negative-index resonance of the structure (see 
Fig. 4 later on). Fabrication tolerances, passivation, and 
quenching may lead to a reduction of the local gain close 
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Figure 2. (color online) Conversion efficiencies over pump du- 
ration Tp for gap sizes of — 10 nm (thick lines: 5 nm). (a) Left 
axis: Total absorption efficiency a (red lines), external pump 
efficiency ar]g (blue lines), and dissipative efficiency arif (green 
lines); absorption coefficient A (yellow squares) obtained from 
ATR calculations; Right axis: Total inversion {N2 — Ni) (black 
lines) of the gain material, (b) Pump intensity Ip over Tp 
indicating the switch-off times (dashed lines), (c) Spatial in- 
version profiles, [N2{r) — A^i(r)] /N, for four different pump 
times (i)-(iv): (top row) x-y cross section in the middle be- 
tween the fishnet metal films; (bottom row) x-z cross-section 
through the centers of the holes. 

to metal surfaces, which we qualitatively accommodate 
for by considering a gap between the metal and the gain 
material with values of g = — 10 nm. The volume V 
used for the rate retrieval [Eq. (2)] is chosen to tightly 
enclose the structure. 

The amplification performance of the active double- 
fishnet critically depends on where and how efficiently 
energy is deposited in the gain medium and the metal 
during the pump process. While macroscopically one 
can only measure the total absorption, the presented 
methodology allows us to trace the conversion of en- 
ergy into inversion and thermal energy microscopically 
within the gain medium and the metal films, as shown 
in Fig. 2. We pump the metamaterial with a pulse of 
constant peak intensity /p = 3.3MW/cm^ but varying 
duration Tp = — 180 ps (in steps of 1.38 ps) and record 
the gradual change in the efficiency with which energy 
is absorbed by the gain medium and the metal. To this 
end, we define the following photon conversion efficien- 
cies: the total absoption efficiency a = — (V • S) / {IpA) 
is the net percentage of photons absorbed in the volume 
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(influx); the internal pump efliciency T^g = Fg/A (with 
Fg = — Fe — Fa) is the percentage of absorbed photons 
transfering their energy to the gain medium; and the in- 
ternal dissipative 'efficiency' r]f = — Ff/A is the percent- 
age of absorbed photons dissipated in the metal. 

Figure 2(a) shows that just after the ramp- up of the 
pump pulse (2.7 ps) approximately a third of all incident 
photons are absorbed inside the cavity (red lines), with 
the efficiency of the energy transfer to the dye (blue lines) 
dominating over dissipative losses (green lines) by a fac- 
tor of r]g/r]f ^ 1.45. The energy absorbed by the dye 
is converted into inversion, which initially concentrates 
around the metal interfaces where the pump fleld en- 
hancement is strong [see Fig. 2(c) (i)]. As the duration 
of the pump increases, the efficiency with which inversion 
is created in the dye {ar]g) decreases rapidly falling below 
the contribution from dissipative losses [Fig. 2(a)]. This 
drop in pump efficiency is caused by the saturation of 
the gain medium, with areas of high inversion gradually 
spreading out from close to the metal, as can be seen in 
Figs. 2(c) (i) — (iv). Eventually, a global saturation of the 
total inversion [Fig. 2(a), black lines] sets in where half 
the maximum value is already reached within the first 
30 ps. 

We further note that the efficiency with which energy is 
absorbed by the gain medium (blue lines) does not reach 
zero for Tp ^ oo since the broad emission line of the dye 
overlaps with the pump wavelength. As the medium ap- 
proaches saturation the internal quantum yield decreases 
and the electric field energy is increasingly absorbed by 
the dye in absorption-emission cycles without creating 
inversion. 

Having pumped the active metamaterial into inversion, 
we now probe it with a signal pulse of duration Tg = 1 ps 
and peak intensity /g = 1.23kW/cm^ (Fig. 3(b), black 
line) at the emission wavelength of the dye. Note that 
the energy inside the resonator (Fig. 3(b), red line) does 
not instantaneously follow the excitation /g, but exhibits 
a rising fiank during switch-on and an exponential decay 
after the switch-off of the probe pulse. After the tran- 
sient build-up of energy in the resonant mode (Fig. 3(a), 
negative values for the black line) two distinct phases 
of constant rates follow: steady-state continuous excita- 
tion (CE) characterized by Ft = 0; and free decay (FD) 
of the excited mode with a Q-factor Q = cJe/^t- Fig- 
ure 3(a) exemplifies the retrieval of rates for a specific 
pump duration (and inversion) and 5nm gap size. By 
performing corresponding probe experiments for pump 
times Tp = — 180 ps and gap sizes ^ = — 10 nm we 
obtain the results reported in Figs. 3(c) and (d). 

In the CE phase [Fig. 3(c)], the energy inside the driven 
resonator does not change with time and, according to 
Eq. (2), the fiux rate A (red lines) then equals the dif- 
ference between the net gain rate Fg (blue lines) and the 
dissipative- loss rate Ff (green lines), i.e. A = Fg — Ff. 
Figure 3(c) shows that overcoming the amplification 
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Figure 3. (color online) (a) Rate dynamics during probe pro- 
cess for 5nm gap in amplifying regime: Gain rate Fg (blue 
line), dissipative loss rate Ff (green line), outflux rate A (red 
line), and energy decay rate Ft (black line), (b) Probe pulse 
intensity Is (black line) and energy U inside the volume (red 
line) over time, with highlighted regimes of continuous exci- 
tation (CE) and free decay (FD). (c,d) Extracted rates over 
total inversion in the gain medium for gap sizes — 10 nm 
under CE (c) and in FD (d) [colors as in (a)]. 



threshold A = 0, where the gain rate equals the loss 
rate, requires higher values of total inversion for larger 
gap sizes. This is because the field of the resonator mode 
at the probe wavelength E(Ae, r), which weighs the local 
gain rate according to Fg oc J Fg(r) |E(Ae, r)|^ d^ [16], is 
highest close to the metal. While the inversion in these 
regions makes larger contributions to the effective gain, 
its impact on the gain-enhancement is less strong than 
one might expect, making amplification possible even for 
larger gap sizes of up to g ^ 7.5 nm. Energy transfer 
close to the metal interfaces (electron- hole coupling [9]) 
is therefore expected to increase the threshold but does 
not necessarily prevent amplification. 

When the probe is switched off, the rates rebalance 
and enter the FD phase [see Fig. 3(d)]. The energy inside 
the resonator then decays exponentially with a rate Ft = 
A + Ff — Fg > 0. Just as the amplification threshold was 
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Figure 4. (color online) Absorption (blue lines) and trans- 
mission spectra (red lines) for pump times Tp = 2.7 — 180 ps 
(lighter to darker colors, black arrows) for 5 nm gap size. The 
Q-factor calculated from the width of the transmission reso- 
nance is given for Tp = 180 ps . Inset: real (blue lines) and 
imaginary parts (red lines) of the effective refractive index 
retrieved from the complex ATR spectra. 



characterized by A = in the CE regirae, we identify the 
lasing threshold as Ft = in the FD regime. Here, for 
lasing to occur the gain would have to overcome the sum 
of dissipative and radiative loss. 

The constancy of the dissipative loss rate (Ff 
22ps~^) and the radiative damping rate in free decay 
(A ^ 48ps~^) show that neither the overlap of the field 
with the metal nor the coupling to the radiative con- 
tinuum are affected by gap size or inversion. This in- 
dicates that the mode profiles in the two regimes (CE 
and FD) remain constant. Interestingly, the dissipative 
loss rate also remains unaltered during the CE FD 
transition, despite the switch-off of the incident field. In 
effect, a window of amplification exists for the gain rate 
Fg between the point of loss-compensation and the lasing 
threshold, Ff < Fg < Ff + A^^. For the considered ac- 
tive plasmonic metamaterial the radiative flux rate A^^ 
(Fig. 3 (d), red lines) is more than two times larger than 
the dissipative- loss rate Ff (green lines). In practice, the 
regime of amplification might be further limited by gain 
saturation as Fg < F|^*. For sufficiently high gain densi- 
ties the lasing threshold can be crossed. The energy in- 
side the resonator then increases exponentially (Ft < 0, 
lasing instability) until gain depletion sets in, which is 
indeed observed. 

The microscopically calculated cavity Q— factor and 
the amplification threshold can be related to the ATR 
spectra depicted in Fig. 4. We see that for longer pump 
durations (increased inversions) all spectra converge due 
to gain saturation, including those of the extracted ef- 
fective refractive index [20] shown in the inset. No 
further increase (decrease) in the transmission (absorp- 



tion) of the probe pulse through the structure is possi- 
ble. The width of the transmission resonance narrows to 
A/ ^ 4.9 THz at saturation, which is in good agreement 
with the extracted decay rate of Ft ~ 31 ps~^ ~ 27rA/ 
[last point of full black line in Fig. 3(d)]. Further, the ab- 
sorption coefficient A at the probe wavelength becomes 
zero at precisely the same inversion value at which A = 
in the CE regime [full red line in Fig. 3(c)]. At the pump 
wavelength, A remains positive, decreases with higher 
values of inversion and, presented as yellow squares in 
Fig. 2, agrees quantitatively with the microscopically cal- 
culated absorption efficiency a. 

In conclusion, we have introduced a generic Maxwell 
Bloch-based methodology that allowed us to investigate 
the dynamic origins of light amplification, i.e. the trans- 
fer of energy between microscopic gain and loss channels, 
in a fishnet metamaterial. The extracted rates associated 
with the subsystems of dye, metal, and radiative contin- 
uum reveal that radiative damping dominates over the 
dissipative loss. Our results show that steady-state net 
amplification free of lasing instabilities is possible over a 
broad parameter regime and may be achieved even when 
gain passivation close to the metal is considered. 

We gratefully acknowledge support by the Leverhulme 
Trust, the Royal Academy of Engineering and the EP- 
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